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Abstract I outline, from a theoretical and somewhat personal perspective, signifi- 
cant features of Pulsar Wind Nebulae as Cosmic Accelerators. I discuss recent stud- 
ies of Pulsar Wind Nebulae (PWNe). I pay special attention to the recently discov- 
ered gamma ray "flares" in the Crab Nebula's emission, focusing on the possibility, 
raised by the observations, that the accelerating electric field exceeds the magnetic 
field, suggesting that reconnection in the persistent current layer (a "current sheet") 
plays a significant role in the behavior of this well studied Pevatron. I address the 
present status of the termination shock model for the particle accelerator that con- 
verts the wind flow energy to the observed non thermal particle spectra, concluding 
that it has a number of major difficulties related to the transverse magnetic geometry 
of the shock wave. I discuss recent work on the inferred pair outflow rates, which 
are in excess of those predicted by existing theories of pair creation, and use those 
results to point out that the consequent mass loading of the wind reduces the wind's 
bulk flow 4-velocity to the point that dissipation of the magnetic field in a pulsar's 
wind upstream of the termination shock is restored to life as a viable model for the 
solution of the "a" problem. I discuss some suggestions that current starvation in 
the current flow supporting the structured ("striped") upstream magnetic field, per- 
haps inducing a transition to superluminal wave propagation. I show that current 
starvation probably does not occur, because those currents are carried in the current 
sheet separating there stripes rather than in the stripes themselves PJ 
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2 Jonathan Arons 

1 Gamma Ray Flares from the Crab Nebula: Observations 



Rotation Powered Pulsars (RPPs) and their Pulsar Wind Nebulae (PWNe) are cos- 
mic accelerators containing the highest energy particles in identified sources. Mod- 
eling photon emission, which is synchrotron radiation from radio (/ieV) up to GeV 
gamma rays, provides evidence for electrons and (probably) positrons with ener- 
gies up to a few PeV, as has been known for many years. In the last two years, the 
new generation of satellite borne gamma ray telescopes, AGILE and FERMI, have 
discovered strong variability in the e > 100 MeV gamma ray flux from the Crab 
Nebula that may revise, or at least extend, ideas about the acceleration of such very 
high energy particles in well-identified sources - the sources of the ultra-high en- 
ergy cosmic rays, whose observed flux contains particles with energy up to ~ 0.1 
ZeV, have yet to be identified (Kotera & Olinto 2010). Figure [T] shows the average 
spectrum of the Crab seen by the Large Area Telescope (LAT) in the first 33 months 
after the launch of the Fermi satellite. 
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Fig. 1 Gamma ray (10 TeV > e > 1 MeV) average spectrum of the Crab Nebula in 2008-1 1 , from 
\Buehler et aL\0T2) . 



The standard model for the high energy emission employs some form of accel- 
eration of non-thermal distributions of electrons (and positrons) at the termination 
"shock" of the relativistic, magnetized radiatively inefficient ("cold") wind that car- 
ries the central pulsar's rotational energy loss from the magnetosphere (dimension 
~ Rl = c/Q* = 1576(P/33 msec) km (the light cylinder distance) to the location 
of the "termination working surfact^j' (the TWS) at Rtws = {Er/ AncP ne b u ia) 1 ^ 2 = 
3.6 x 10 12 km =0.4 Light years = 0.12 pc = 2.3 x \Q 9 R L . For the Crab Pulsar, tim- 
ing observations of the pulsar" period and period derivative yield Er = I£2*Q t = 

1 I prefer this somewhat old fashioned term for the surface or region where outflow energy converts 
to non-thermal heat, since as will be seen, the interpretation of the surface as a traditional MHD 
shock wave does not work, requiring one to endow the word "shock" with unusual properties, if it 
is to be used in this context. 
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5 x 10 38 ergs/s, with a factor ~ 2 uncertainty due to uncertainty in the magnitude 
of the neutron star's moment of inertia /, itself due to imperfect knowledge of the 
equation of state of the nuclear matter in the star's interior. 

Figure [2] shows the Nebula's light curve above 100 MeV using 12 hour aver- 
ages during the same period. Apparent are the restless small amplitude variability 
and the three flares (February 2009, September 2010 and April 201 1) seen by AG- 
ILE and FERMI. Figure [5] shows the same kind of data for the most spectacular 
event observed in April 201 1. During that event, the peak power reached 1% of the 
pulsar's spin down power (the ultimate source of the energy), thirty times the aver- 
age gamma ray brightness at the same energy, with factors of a few flux variation 
on time scales of hours (Bal bo et al.\ |20lT] l. The 2011 flare spectrum drastically 
hardened compared to the normal, average nebular gamma ray spectrum, as seen in 
Figure [4] Despite the restless behavior of the y-ray emission at all times, the events 
identified as flares do seem to be special events - unlike the largest earthquake^] 
they are not the tail of a (non-Gaussian) distribution of amplitudes (J. Scargle, 201 1, 
personal communication). Throughout all these events, the pulsed emission coming 
directly from the pulsar did not change - no variations in the pulsar spin aside from 
the regular spin down have been detected, and no changes in the spectra or light 
curves of the pulsed emission appeared (Abdo et al. 2011 Buehler et al. 2012[ ). 



Thus, the blame for these phenomena has to be assigned to events in the nebula, 
the calorimeter which makes luminous the otherwise dark outflow from the central 
compact object. 

The wind model for the excitation of the Crab Nebula, introduced in a form 



mixed with the vacuum wave model for pulsar spin down by Rees & Gunn ( 1974 1 



and re-introduced in its fully MHD form by Kennel & Coroniti] ( |1984]l^ is designed 
to quantity the steady emission from the Crab Nebula and other PWNq^] It has been 
shown to incorporate unsteady, quasi-periodic variability with time scale ~ a year 
(Spitkovsky & Arons 2004 Camus et al. 2009) ~ flow transit time from pulsar 
to the TWS and in the post TWS flow, with simulated behavior rather akin to the 
variable "wisp" features long known fro m optical (|Lampland|[T92T[|Scargle 1969 1, 
X-ray fWeisskopf et qZ.| |2000 ) and radio ( jBietenholz et aZ.||2001||2004| l studies. The 
wisp variability time scale's similarity to the time separating the well demarcated 
flares perhaps suggests the origin of the gamma ray flares lies in the physics of 
the TWS. As it stands, however, that model has no ready explanation for the short 
variability time during a flare - the hierarchy of time scales (many months between 
flares, few days flare duration, and few hours intra-flare variability) has not been 
part of the modeling to date. In addition, the spectral hardening, looking as if a new, 
transient component of the PeV particle distributions appeared, then vanished, poses 
a serious challenge to the standard model. I discuss the significance of these results 
in m 



Carlson, Langer & Shaw ( 1994) for a review of earthquake dynamics. 



' The relativistic MHD wind model for pulsar spin down was introduced by Michel > 1969 
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Fig. 2 Gamma ray (e > 100 MeV) light curve of the Crab Nebula as detected by the Fermi LAT, 
for the thirty-three months ending in April 201 1, binned in 12 hour intervals (Buehler etal. 2012 ). 
The flare of September-October 2007, detected by AGILE before Fermi's launch l |Tavani it et al.| 
201 1 ), is not shown here. 
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Fig. 3 Gamma ray (e > 100 MeV) April 201 1 flare light curve, from 
the fast (few hours) variations . 
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Fig. 4 Gamma ray (e > 100 MeV) April 2011 flare spectra at 11 times during the flare, from 
Buehler et al. {2012 1. The time windows are indicated by the numbers in the bottom left corner 
of each panel and correspond to the intervals shown in Figure '3 The dotted line shows the SED 
of the flaring component (logarithmic flux eF E in ergs cm~ 2 sec -1 as a function of photon energy 
e), the dot-dashed line the constant background from the synchrotron nebula, and the dashed line 
is the sum of both components. The average Crab nebular spectrum in the first 33 months of 
Fermi observations is also shown in gray for comparison. Superficially, these spectra suggest the 
acceleration of an almost monoenergetic component of the radiating particles (a "beam"?), with 
this unusual accelerator persisting for a few days, or, if lasting longer, beamed toward us for only 
a few days. 



2 Pulsar Wind Nebulae: Basic Concepts and the Standard Model 

PWNe and their generative pulsars provide the first known and most definitive ex- 
ample of compact astrophysical systems which draw the power for their observed 
emissions from the electromagnetic extraction of rotational energy from gravitation- 
ally bound objects. They have motivated models for similar energy extraction from 
disks around other gravitating bodies, such as black holes (e.g. Rees|198 4 Begel- 



man et al. T984J). As objects of study, the RPPs have a virtue lacking in other systems 



studied in high energy astrophysics: timing of the precisely measured pulse periods, 
uniquely interpretable as the rotation periods of the underlying neutron stars, pro- 
vide measurements of the total energy budget free of all astrophysical ambiguities, 
other than the factor of ~ 2 uncertainty in neutron stars' moments of inertia, arising 
from the lack of precise knowledge in the equation of state of dense matter. The 
measured rate of rotational energy loss, 
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E R = -mA = 4nIp, (1) 

tells us the total energy budget for these systems, without our having to understand 
anything about the photon emissions. The observed distributions of P,P appea in 
Figure[5] This is both a blessing and a curse - a blessing, because in contrast to other 
relativistic astrophysical systems, we know the total luminosity in all forms, seen 
and unseen, without having to unravel the partition between flow kinetic energy, 
large scale Poynting flux, thermal energy and radiative losses - a curse, because the 
energy loss is radiatively silent, thus supplying little information as to the details of 
the energy outflow, leaving the mechanics of the machine mysterious. 



2.1 Magnetospheres 



Nevertheless, progress has been made on interpreting the physics of the observed 
radiation. The most significant advances have come in interpreting the Pulsar Wind 
Nebulae (PWNe), since these act as catch basins for the rotational energy lost. Ob- 
servations of these systems have made clear that RPPs deliver their energy to the 
outside world in the form of highly relativistic, magnetized outflows - stellar winds 
that are exaggerated versions of the solar wind - which must be electromagnetically 
driven by the magnetic pressure of the wound up magnetic field. The strength of that 
field is estimated by using the theory of magnetic braking of the neutron stars' spin, 
which suggests 

G = -KG", G =2n/P, (2) 

applied to the observed rotation periods and spindown rates P = —2nG/G 2 . 

The earliest model applied vacuum electrodynamics to a rotating sphere endowed 
with a magnetic dipole moment fl centered at the stars' centers and tipped with 
respect to the rotation axis by an angle i. That theory yields the spindown lumi - 



1967 



Ostriker & Gunn 1969). 



nosity E R = KG 4 , K = (2/3)^ 2 sin z //c 3 ( |Pacini 
Thus n = 3 in this model, if fl and / are constant. Vacuum theory was motivated 
by the large gravitational forces at the surfaces, suggesting no plasma more than a 

( Hoyle et al.\ 1964). But immediately the much earlier 



meter or so above the star 



observation (Deutsch 1955 1, made in the context of magnetic A stars, that in vac- 
uum large electric fields parallel to B would overwhelm gravity and pull charged 
particles out from the star until the vacuum electric field would be altered, reduc- 
ing E ■ B down to zero, was recovered ( |Goldreich & Julian] |1969| l and extended 
with the suggestion that the charged particles would feed a curious charge sepa- 
rated wind, with a total electric current / = c<J> mag , ^mag = total magnetospheric 
potential = y/E R /c = 4 x 10 16 (P/3 x 10-13) 1 / 2 (33 msec/P 3 ) 3 / 2 Volts - the nu- 
merical result is normalized to the Crab pulsar. That wind could carry away the 
rotational energy of even the aligned rotator, in a Poynting flux dominated flow 
- the electromagnetic energy density would vastly exceed the kinetic energy den- 
sity (and pressure) of the outflow, in the initially conceived model - the particle 
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flux in that scheme is only c<P/e = 2.3 x 10 30 (/45f > i5/^ >3 ) 1 ^ 2 elementary charges/s, 
Z45 = 7/10 45 cgs.As = P/10~ 15 , the "Goldreich-Julian" current. 

The charge separated model has several really serious theoretical difficulties, but 
perhaps of greater importance is that the observations of young PWNe (<£ ma g > 
10 14 5 Volts) have made clear that the particle outflows are many orders of magni- 
tude larger than the Goldreich-Julian current, thus motivating MHD models, where 
E ■ B = is assumed from the start. Pair creation occurring somewhere within the 
magnetospheres (first suggested by Sturrock 1971 1 is thought to be the origin of 



the required dense plasma, although a quantitative model that yields the high mass 
loss rates observed is still lacking (Bucciantini et al. 20 1 1 1 and references therein.) 
In recent years, solutions for the MHD structure of the magnetosphere in the ap- 
propriate force-free limit have been obtained (numerically - attempts of analytically 
minded theorists to guess the detailed answer uniformly failed over 30 years of try- 



ing), first for the aligned rotator ( Contopoulos et al. 1999 Gruzinov 2005 



sarov 



Komis- 



[2006} |Timokhin| |20"06), then f or the full 3D oblique rotator (Spitkovsky 



2006 Kalapotharakos & Contopoulos, 2009). 



As far as the P, P diagram goes, the main inference from modeling the force-free 
magnetosphere is the innocuous looking result for the scale factor K in expression 



3, ,2 



K = k{\ 



sin 1) — 5- 



k= 1 ±0.1. 



(3) 



Physically, the most important result has been the identification of the current sheets 
separating the closed from the open regions of the magnetosphere, extending into the 
wind beyond the magnetosphere, whose last closed flux surface ends just touching 
the light cylinder whose cylindrical radius is CJ = c/Q. That such a current sheet 
should be present has been suspected from the early days of RPP research (e.g., 
Michel 1975). The error in k reflects the uncertainties in the numerical treatment 
of the problem, many of which are associated with how the current sheet is repre- 
sented in the numerical schemes. Figure [6] shows a slice through the 3d force-free 
magnetosphere of the 60° rotator, from Spitkovsky's (2006) results. 
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Fig. 5 Observed RPP periods and period derivatives, adapted from Kaspi 1 2010 . The line bound- 
ing the pulsar population corresponds to the magnetospheric voltage 4> mag = 10 12 V. It clearly 
marks a boundary beyond which pulsar emission is unlikely. Simple estimates of pair creation 
within the magnetosphere, which assume the existence of a parallel voltage drop in the strong field 
region r < Rl with magnitude a large fraction (generally greater than a few %) of <£ mag , suggests 
this source of plasma should occur only for <P mag > 10 12 V. The numerical value is determined by 
the pair conversion opacity of photons propagating in the super strong B field and the magnetic 
geometry; see Hibsch man & Arons| ( |200l) and Medin & Lai (2010) for modern evaluations. This 
result underpins the idea that pair creation is essential for radio emission jSturrock] |1971} . The 
change in slope of this "death line" at short period, small P indicates more sophistication in the 
pair creation physics, and/or in the association of pairs with radio emission, than is incorporated in 
the simplest models, a conclusion also apparent from the variety of quantitative problems with this 
widely accepted hypothesis (Hibschman & Arons, 2001 Medin & Lai, 2010 1. 



As far as the P, P diagram goes, the main inference from modeling the force-free 
magnetosphere is the innocuous looking result for the scale factor K in expression 
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K = k(l+sm 2 i) — f-,/t=l±0.1. (4) 

Physically, the most important result has been the identification of the current sheets 
separating the closed from the open regions of the magnetosphere, extending into the 
wind beyond the magnetosphere, whose last closed flux surface ends just touching 
the light cylinder whose cylindrical radius is GJ = c/Q. That such a current sheet 
should be present has been suspected from the early days of RPP research (e.g., 
Michel 1975 ). The error in k reflects the uncertainties in the numerical treatment 
of the problem, many of which are associated with how the current sheet is repre- 
sented in the numerical schemes. Figure [6] shows a slice through the 3d force-free 
magnetosphere of the 60° rotator, from Spitkovsky's (2006) results. 




olar Gap 



'Slot Gap 



Outer Gap 



Fig. 6 Field lines and current density of the oblique force-free rotator, ( = 60°, from Spitkovsky 
(2006), seen in a cut in the il , fl plane. The last closed field lines end at a Y-point (a Y-line, in 3D) at 
distance Ry from the neutron star - technically, the closed zone ends at cusp, not a Y-line, since the 
separatrix between the closed and open zones interior to the line of closure carries the return current 
Uzdensky & Kulsrud 1997 but this distinction is of no importance to the present discussion. The 
current sheet encloses and separates the closed from the open field regions of the magnetosphere, 
and the separate branches merge in the wind zone, where the folded sheet continues to separate 
the oppositely directed fields of the striped wind. The arrows and "Gap" labels locate sites where 
vacuum gaps have been postulated, in test particle models of accelerators that lead to gamma ray 
emission and pair creation. None of these gaps appear in a current sheet accelerator based on the 
force-free magnetosphere model. 
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Observationally, expression Q shows that the vacuum rotator's "braking index" 
— 3 is preserved in full force-free MHD, which contradicts the observed values in 



the small number of stars where n has been determined ( |Livingstone et al. 2007 ), 
a contradiction which has led to a variety of suggestions ranging from evolution of 
the magnetic moment ji or the obliquity / (Blandford & Romani 1988)1 to effects of 



reconnection on the rate of conversion of open magnetic flux to closed ( Contopoulos| 
& Spitkovsky 2006), as the star spins down and the closed zone expands at the 
expense of the amount of open magnetic flux. That reconnection might affect the 
braking index is readily derived from the fact that the torque really depends on the 
magnitude of the open magnetic flux.The amount of open flux depends on the size 
of the closed zone, which ends at Ry. If Ry /Rl < 1, the torque increases because of 
more open field lines and larger Poynting flux than is the case for a magnetosphere 
closing at r = R^. Bucciantini et al. (2006) show that the braking index is 



QQ 



din 



:3+2- 



(l+§) 



<?lni2 



(5) 



thus, If Ry/Rt decreases with decreasing £2, then n < 3. Reconnection usually is 
unsteady - Figure [7] shows the blobs ("plasmoids") ejected from the Y-line at the 
base of the current sheet in the relativistic wind (beyond the light cylinder) of the 
aligned rotator. 




Fig. 7 Plasmoids formed at the base of the current sheet of the relativistic aligned rotator, where 
the sheet crosses the light cylinder, from a relativistic MHD simulation of a newly born neutron 
star's magnetosphere (Bucciantini et al. . 2006). The dissipation that allows reconnection to occur 
is numerical. They move out radially at the local Alfven speed sa c, and recur on the magneto- 
spheric Alfven transit time ~ P jiz 
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The speculation is that spindown slightly biases these reconnection events so that 
on the much longer spindown timescale, the net open flux slowly converts to closed. 
These current fluctuations might be associated with the timing "noise" (|Arons 



1981a Cheng 1987|l iden tified long ago with torque fluctuations (e.g., JHelfand et 
al. 1 1980 Scott et a/.|2003 1, although recent analysis of longer data sets ( |Lyne et al 



2010|l in long period pulsars has called the pure noise interpretation into question. 



In itself, the force-free model does not provide mechanisms for photon emission. 
But it has a variety of implications, which are slowly being addressed. 

• The model specifies the polar flux tube size and shape - for i ^ 0, it is noncircular 
with a polar cap center displaced from the magnetic axis, even when the mag- 
netic field is the simplest, that of a star centered, point dipole (Bai & Spit kovsky] 



2010). This has consequences for radio polarization structure, and for polar cap 



areas and dipole offsets inferred from soft X-ray emission from polar caps (Bog- 
danov et al. |2007[), thought to be heated by magnetospheric particle bombard- 



ment (Arons M98Tb Zavlin & Pavlov 1998 Harding & Muslimov 2002). These 



theoretical improvements of the polar cap model have yet to be noticed and incor- 
porated in phenomenological models of the observations used by data analysts; 
such incorporation might yield interesting tests of the force-free model. |Chung| 
& Melato"sl ( |2011| ) show how "Stokes Tomography" applied to radio polarization 



might be usefully employed in this task. 
The force free model quantitatively specifies the return currents required to pre- 
vent the star from charging up, as the polar flow extracts charge from the star. 
The results for the oblique rotator are partly in accord with long held expecta- 
tions, that in part return current exists in a thin ("auroral") sheet bounding the 
polar flux tube (e.g. Goldreich & Julian 1 969] |Michel| 1 975" i, consistent with the 
open circuited model (Goldreich & Julian| 1 1969 ) - current closure occurs far 



away, in the nebula/interstellar medium beyond the wind termination shock or 
perhaps in the outer wind - plus qualitatively new features: a) part of the return 
current surrounding the polar flux tube is spatially distributed, even in the aligned 
rotator; b) in the oblique rotator, part of this current system is not a return cur- 



rent at all, but couples the two polar regions together (Bai & Spitkovsky 2010 1 
- in the orthogonal rotator (i = 90°), the auroral component of the current is en- 
tirely in the polar coupling flow, with the volume current out of each half of the 
polar cap having equal amount and opposite sign, also consistent with early ex- 
pectations (e.g. Scharlem ann et. a/.|1978| ) - the orthogonal rotator automatically 
balances its charge loss. A preliminary account of some aspects of "auroral" par- 
ticle acceleration based on the force-free model are outlined in Arons (201 1 ). The 
radiative consequences of these features are as yet unexplored - for example, the 
spatially distributed part of the return current might be a good candidate for the 
site of "conal" component of pulsar radio emission, an idea which requires non- 
force-free modeling of the current flow and identification of a workable emission 
process within that current flow modej^Jbefore one can relate the theoretical force 



4 Some earlier ideas on this subject relating to field aligned acceleration and gamma ray emission 
can be found in |Arons||1983b| , |Gruzinov||2008} , for example. 
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free current distributions to the observations in a testable manner (although easier 
kinematic comparisons are certainly possible). 

The location of the return current layer having been determined, the hypothesis 
that the return current layer is the site of the beamed particle accelerator that 
gives rise to the pulsed gamma rays observed by the FERMI and earlier orbiting 
gamma ray telescopes (see |Ray & Saz Parkinson|201 1} can now be investigated 
in the context of a self consistent magnetospheric structure that allows a quanti- 
tative evaluation of the beaming characteristics implied by the radiating current 
sheet concept - see Bai & Spitkovsky (2010) for a kinematical study of the radi- 



ating current sheet idea. 



2.2 Current Sheet and Pulsed y-ray emission in the Lighthouse 
Model: Relativistic Aurorae 

Making progress on a physical model for radiation from the current layers can most 
expeditiously take advantage of the facts that a) pulsed gamma ray emission, when 
observed, is the largest photon output from rotation powered pulsars, but b) gen- 
erally has less luminosity than the spin-down luminosity of these stars. Figure [8] 
illustrates this fact, which summarizes the results from the LAT instrument as of 
Spring 2010. Thus the energy invested in particle acceleration is a small fraction of 
the energy stored in Poynting fluxes, for radiation reaction limited acceleration, so 
that the force free model can be considered as a good zeroth order magnetospheric 
description. 

The voltage limitation prediction L gamma °< ^mag W is the same as Ly °< "parti- 
cle current" = particle flux = Goldreich-Julian flux in unidirectional beam models of 
the polar electric current flow (Harding 1981) only if the accelerator carries a fixed 
fraction of the total electric current that enters into the spin-down torque, indepen- 
dent of P,P, In the traditional slot or outer gap models, pair creation establishes the 



limitation of the accelerator to a thin sheet either in the outer magnetosphere (Cheng, 



Ho & Ruderman 1986) or back in the polar cap (Arons 1983a I. In such models, 



the gap width w* of a model that successfully reproduces the sharply peaked light 
curves, projected onto the neutron star following the field lines, is necessarily small 
compared to the polar cap size, and varies with the pulsars' magnetic moment and 
spin parameters. Thus Ly °< <£ ma g is not the same as Ly °< electric current / = c<P ma o. 
This is not an issue for a model based on the currents flowing in the return current 
layer, which necessarily carries the whole magnetospheric return current, for most 
obliquities of the magnetic moment with respect to the rotation axis. 

That magnetospheric current sheets, with particle densities within the sheets high 
compared to the Goldreich-Julian value, can sustain large parallel electric fields is 
well known in planetary magnetospheres - such sheets are the accelerator sites of the 
particle beams that stimulate the aurora observed in the upper atmospheres of the 



Pulsar Wind Nebulae as Cosmic Pevatrons: A Current Sheet's Tale* 



13 



Earth, Jupiter and Saturn, for exampl^An elementary illustration of this possibility 
comes from considering the inertial term in the generalized Ohm's for the electric 
field parallel to B, which can be written as, in the relativistic case, as 

4n ( dJn B „ . ,, ,..1 my I 

•V-[y(Jv+vJ)Uoc-L_ , (6) 



" 0)j \ dt B ' J n 4 C urrentPfl 

with my the particles' relativistic mass, n their density, / the total current set by the 
force free magnetosphere, ^current is the thickness of the current carrying channel, 
and Pb is the radius of curvature of the magnetic field. 
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Fig. 8 Ratio of observed gamma ray luminosity Ly to observed spindown power Er = L S pindown for 
LAT pulsars ( Abdo et al. 2010 ). If the acceleration of the gamma ray emitting particles is radiation 
reaction limited (generally true in models that assume curvature radiation as the emission mecha- 
nism, as is the case in most of the "gap" models), Ly is a good proxy for the acceleration power. The 

■ —1/2 

trend at higher luminosities roughly follows the efficiency Ly/ Z. S pindown x E R expected if some 
mechanism limits the accelerating voltage A <P to a fixed fract ion of the magnetospheric voltage 
4> mag = \JHgJc, (e.g., |Arons| 1996 Harding & Muslimov|2002| , a limitation usually (and plausi- 
bly) attributed to pair creation. The large dispersion in Figure[8]comes primarily from the uncertain 
distances, and also from the uncertain beaming correction required to go from the fraction of the 
sky illuminated by the beam to the total radiative - the LAT team assumed uniform phase-averaged 
beaming across the sky (1 sterradian). 



Thus inertia of the current carriers can act as a effective resistance in these high in- 
ductance systems, which establish the currents electromagnetically. The presence of 
this "resistor" in the circuit forces a parallel electric field to appear. The same is true 
of any parallel (to B) load - pressure and radiation reaction are other effects which 
can serve as loads, with pressure especially important in the diffusion region around 



For a review of such phenomena, see 



Paschmann et al. 



2002 i. Of particular significance to 



pulsars is the fact that the field aligned currents that power narrow auroral arcs consist of pre- 
cipitating electron beams launched from the reconnection region in the distant magnetotail and 
counterstreaming ions launched from the planetary atmosphere. 



14 



Jonathan Arons 



the singular line shown in figure [9] Equation |6]) is most useful in the corotating 
frame, and when the current is due to relative motion between the species (elec- 
trons and positrons, and in some circumstances heavy ions) that is slow compared 
to the bulk fluid velocity v. In the current circumstance, it turns out that the current 
is better described as counter-streaming beams - in that case, describing the beams 
as separate fluids is more appropriate, and they can form the total plasma density 
in the current flow channel, rather than being a low density component in a much 
denser plasma. 

Expression (|6| does make clear that acceleration is prone to maximize when the 
the relativistic mass is high and the current density is high (large //^current) - ^current 
is generally microscopic, expected to be on the order of c/(O p , and is established 
by the dynamics of the the singular region where the closed zone ends, illustrated 
(in the cartoon approximation) in Figure [9] The capture rate of pair plasma into the 
diffusion region is 



2l D A L /3 rec c<£ mag 

— 2 



Vrm ^ rrec —mag .„ 

iv diflusion± D 2 o ± - ' ' ' 



with K± the multiplicity (multiplier of a fiducial Goldreich- Julian outflow rate 
c<5 ma g/e) that gives the number of pairs in the total outflow, jSrec the reconnection 
speed in units of c and j3 w ; n( i the polar wind outflow velocity (set equal to unity). Al 
is the thickness of the current channel at the light cylinder, assumed equal to the half 
height of the diffusion region, and Id is the length of the diffusion region. Pressure 
in the diffusion region expels the captured pairs from the diffusion region, outwards 
along the current sheet in the wind and inwards along the auroral current channels. 
Expression (|7jl assumes the plasma flux in the wind has a gradient across B, as is 
likely since the accelerating electric field in the polar cap that leads to the pair cre- 
ation that feeds the wind is small near the cap edge. The pressure supported electric 
field provides the accelerator which sorts the particles in the diffusion region into a 
precipitating beam (electrons in the geometry shown in Figure |9]l and an oppositely 
charged beam traveling outwards in the wind current sheet - the charge signs of the 
beams are as required by the global electrodynamics. The channel thickness is al- 
most certainly comparable to the skin depth in the pairs. Taking the plasma gradient 
into account leads to the lower limit to the channel thickness 

= rg \ , c „ « R L ( V3 , (8) 



(D p ±[8 = c/co p ±(8)} \2K ± e<Z> x 

where 8 is the distance across B from the formal current sheet location, with a 
precipitating beam flux in the return current channel 

Id c& mag /e f R L X ' 



y± is a measure of the flow momentum and of the comparable momentum dispersion 
of the polar plasma flow emerging from the inner magnetosphere, predicted by pair 
cascade models to be on the order of 10 2 . Numerically, (18JI yields a very small value, 
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on the order of meters to hundreds of meters, the specific value depending on ^mag 
and K±. The kinetics of relativistic reconnection being a largely untrodden subject, 
Id is more or less unknown - it could be as small as 4 £ itself (Petscheck style recon- 
nection), or as large as appears in the numerical dissipation driven reconnection ob- 
served in the force free simulations, Id ~ O.IRl (A. Spitkovsky, personal communi- 
cation). In the terrestrial magnetosphere, satellite observations suggest the diffusion 
region length is intermediate between the ion skin depth and the macroscopic scales. 
For the discussion here (and in the more detailed report in preparation), treating Id 
as a parameter to be constrained by model comparisons to observations appears to 
be the wisest strategy. In analogy to observed non-relativistic reconnection, one ex- 
pects j3 lec ~ O.lvAifven/c = 0.1, a value supported by the few PIC simulations of 
relativistic reconnection [e.g., |Zenitani & Ho shino ( 2 007) >] . 



The precipitating particles form the "hanging charge clouds" invoked by Goldre- 
|ich & Jurian| (fl969 ) to cause the electrostatic extraction of return current from the 
star's atmosphere, which happens if riy(R*) = Fy(R*)/c ^> ncj(R*), or , from 
if 

Id Ae 



JC±»1. (10) 



'wind 



The condition in ( 10 1 is satisfied if Id is large compared to the skin depth, but still 
small compared to the numerical dissipation determined length observed in the force 
free simulations. It can be shown ( Arons|201 1 and in preparation) that the conse- 



quent minimum total potential drop in the twin beam channel is 

1 /3 

. , 1 , R* ( m ± c 2 Y± Avind\ h 

A&mm » -r*m aff — „ — COSi 



mag tf L \2e<P mag K ± ftec J A L 

1 ^ R* ( Avind \ l ^ h 
Rl \ A-ec / Rl 



o^magTr I ~5 I -^COSi. (11) 



If Id/Rl really is as large as 0.1, rather than comparable to the minimum scale Ai, 
the accelerating potential in the return current channel is more than large enough 
(around 10 13 Volts, in the Crab and Vela pulsars) to drive curvature gamma ray 
emission from the beams, which limits the particle energies by radiation reaction. 
Such energies are high enough to lead to pair creation, which may limit the acceler- 
ation, although as is clear from expression parallel potential drops in the current 
carrying region can be sustained even if the plasma is dense. 

The high energy radiation might be synchrotron emission. The counterstreaming 
beams are electromagnetically two stream and shear unstable. Since (O p b eam can 
be comparable to the relativistic cyclotron frequency in the outer magnetosphere, 
the growing waves can excite finite Larmor gyration of the particles in the current 
channel, thus producing incoherent emission through hard X-rays and gamma rays 
are possible. These X-rays are an alternative to soft photons from the star, as targets 
for 7— 7 pair production. If the lower frequency waves can escape the plasma, they 
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Fig. 9 Upper Panel: Electric current structure of the oblique force-free magnetosphere, inclination 
= 60°, from Spitkovsky (2006 1. The current sheet, indicated by the darker color, bounds the closed 
zone. The closed zone ends at a singular Y-line (in the ideal force-free electrodynamics approxi- 
mation) at the light cylinder distance from the neutron star. Lower Panel, left: Current structure at 
the polar cap, illustrated for acute angle between the magnetic moment and the angular velocity 
(electron polar current, £2 ■ jU > 0). The return current in the current sheet consists of a precipitat- 
ing electron beam, launched from the diffusion region around the Y-line and possibly augmented 
by high altitude pair creation within the current sheet, with the charges in the precipitating beam 
extracted by reconnection flow from the pair plasma flowing from the polar cap into the wind, plus 
a counterstreaming ion beam extracted electrostatically from the stellar atmosphere. In obtuse ge- 
ometry (ion polar current, £1 ■ JU < 0), positrons precipitate from the Y-line and counterstreaming 
electrons are extracted from the atmosphere. For clarity, the part of the return current not contained 
in the current sheet is omitted, even though this part of the current system is of increasing signifi- 
cance as i — > 90°. Lower Panel, right: Possible structure of the Y-line region, with the termination 
of the closed zone to the left and the merger of the winds from the two polar caps to the right. 
The "guide field" 6^ also reverses across the mid-plane of the flow, along with the poloidal open 
field. As reconnection occurs, some of the ouflowing plasma (speed c/3, /3 rs 1) deflects toward the 
singular, unmagnetized "diffusion" region around the Y-line with speed v rec ~ O.li'Aifven = 0.1c. 
The figure represents a steady (in the co-rotating frame) flow model - in reality, the reconnection 
is likely to be bursty, as in Figure[7] with formation of sporadic X-lines. 
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2.3 Follow the Mass 



The Force-Free model says nothing about the nature of the plasma, that localizes 
£j[ to a boundary layer in the magnetosphere and in the wind that emerges. Pul- 
sar Wind Nebulae (PWNe) demonstrate that pulsars lose rest mass at a rate large 
compared to the fiducial electrodynamic particle loss rate c<P maz /e [for a review 



of PWNe phenomenology, see Gaensler & Slane (2006 1]. The only known expla- 
nation is pair creation in the pulsars' magnetospheres. In many nebulae, the X-ray 
emitting particles rapidly lose energy to synchrotron radiation. Then the nebulae are 
particle and energy calorimeters, allowing direct inference of the pair multiplicity 
in the wind, of TeV to PeV pairs. The measured injection rates of X- and y-ray 
emitting particles N±, up to ~ 10 385 pairs/s, compare well to the predictions of ex- 
isting pair creation models (e.g. Hibschman & Arons|2001| l, yielding multiplicities 
K± = eN±/c<P mag up to ~ 10 4 . However, PWNe are also radio synchrotron emit- 
ters, radiation that samples much lower energy populations (100 MeV to 10 GeV), 
whose radiative efficiency is much less than their X-ray emitting cousins. The re- 
sult is a much larger population of pairs, whose radiative lifetime exceeds that of 
the nebulae. The most efficient hypothesis is that these particles come from the em- 
bedded pulsars also, an idea supported by spectral continuity and by exotica, such 
as the observation of radio "wisps" near the Crab pulsars (Bietenho lz et al.\\200\\ . 
Applying simple evolutionary models allows one to infer time averaged injection 
rates. Early estimates Shklovsky ( 1968) and more recent evaluations by e.g. Atoyan 



& Aharonian] (l996 ) 



Slane et al. (20101, Bucciantini et al. ( 201 1] > yield lower Urn 



its for multiplicities K± all in excess of 10 3 and upper limits for wind 4-velocities 
Jwind = Er/Mc 2 = e<P mag /2K±m±c 2 all less than 10 5 in a number of nebulae. The 



data are the best for the younger systems, although even for these, the lack of far 
infrared data inhibits the analysis. 



PWN Name 


<£mag (PV) 


Age (yr) 




-^wind 


Crab 


100 


955 


10 6 


5 x 10 4 


3C58 


15 


2100 


io 4 - 7 


3 x 10 4 


B1509 


120 


1570 


10 5.3 


1 x 10 4 


Kes 75 


22 


650 


10 5 


7 x 10 4 


W44 


1 


20.3 x 10 3 


10 5 


10 4 


K2/3 Kookaburra 


5.5 


13 x 10 3 


10 5 


10 4 


HESS J1640-465 


3.5 


10 4 


10 6 


1Q 3.6 



The inferred multiplicity excesses are a puzzle for pair creation theory, which 
have been apparent since Shklovsky| ([1968 ) inferred a total injection rate of 10 41 
electrons/s needed to supply the radio emission of the Crab Nebula. These may be 
resolvable by appealing to magnetic anomalies near the neutron stars' surfaces, the 
simplest being an offset of the dipole center from the stellar center, which strength- 
ens the magnetic field at one pole (Arons 1998 Harding & Muslimov||20ir a&b). 
The increase this gives to the magnetic opacity can be greatly enhanced if the mag- 
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netic axis is also tipped with respect to the radial direction, since then gravitational 
bending of photon orbits with respect to the B field direction much increases the 
magnetic opacity for pair creation. Such phenomenological modifications of the low 
altitude magnetic field must respect the observation that radio beaming morphology 
is consistent with the magnetic field being that of a star centered dipole quite close 
to the star ( |Rankin||l990| [Kramer et al]|1998| ). This problem warrants quantitative 
investigation. 

The large inferred multiplicities imply the wind 4-velocities r^j na - to be small 



compared to the much quoted value of 10 6 inferred by Kennel & Coroniti ( 1984 1 in 
their model of the Crab Nebula's optical and harder emission. The large mass load- 
ing and the inferred low wind four velocity has a large impact on the much storied 
"(J problem" of pulsar winds. In ideal MHD, the ratio a = B 2 / 8nm±n±r w j„dC 2 of 
magnetic energy to kinetic energy in the wind is conserved outside the fast mag- 
netosonic radius (since for a cold flow the wind does not substantially accelerate 
outside this surface) and is large - even with the increased mass loading found from 
recent nebular studies, a is always well in excess of several hundred. Nevertheless 
the wind behaves at its termination shock as if O at that distance is small - MHD 
models of the nebulae suggest a at the termination shock is on the order of 0.02 in 
the Crab Nebula (e.g |DelZanna et a/.|20 04) and similar values are plausible in other 



systems. 

Coroniti (1990) suggested that because the wind of an oblique rotator has the 



magnetically striped structure shown in Figure [TOJ magnetic dissipation of the cor- 
rugated B field, generically of a resistive nature propelled by instabilities of the cur- 
rent flow in the current sheet separating the stripeajmight destroy the magnetic field 
of the wind interior to the termination shock, thus converting a high a flow into a 
weakly unmagnetized plasma (a <C 1). If the current sheets separating the magnetic 
stripes are to merge with a speed v s < c as measured in the proper frame of the flow, 
before they reach the termination shock located at distance Rts from the neutron 
star, the merger time in the PWN frame T meige = nr 2 ind (R^ / v s ) must be less than the 
flow time to the termination shock Rjg/c, therefore r^,i n( j < ^{R ts /kR l )(v s /c) = 
5 x 10 4 -\/v s /c (Crab) must be satisfied if Coroniti's model is to be viable - see 
|Lyubarsky & Krrk|pUoT) , |Kirk & Skjaeraasen| ( [2003l l and | Arons | ( |2008> for models 
with greater or lesser detail that address this issue. This inequality is satisfied for 
multiplicities above 10 s , which does appear to be the case for the young PWNe re- 
cently analyzed Bucciantini et al. ( 201 1) . Possible mechanisms that can lead to the 



necessary dissipation are collisionless tearing and drift-kink instabilities of the cur- 
rent sheet, considered as if it were a flat sheet ( |Coroniti| [1990 Zenit ani & Hos hino , 
2007)1 a nd an interesting Weibel-like instability due to interaction between the sheets 



( Arons 2008 1, an effect strongest in the equatorial sector where the folded sheet ap- 
pears locally as neighboring flat sheets with antiparallel current flow in the latitude 
direction. 



6 This wrinkled current sheet, frozen into the wind, is the continuation of the sheet separating the 
closed and open zones interior to the light cylinder, as is apparent in Figure|9| 
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Such effects can be enhanced if the acceleration in the wind zone Rl< r < Rtv/Sj 
driven by the pressure gradient in the wind, itself created by the dissipative heating 
of the wind (Lyuba rsky & Kirk] [2001), drives Rayleigh-Taylor instabilities of the 
current sheets (Lyubarsky, 2010). If this occurs, the small scale disruptions of the 
current sheets can enhance the rate of magnetic decay, by shrinking the spatial scales 
over which the dissipation has to operate. A probable limitation of this idea is that 
it is formulated in such a way as to neglect the magnetic pressure gradient between 
the center of each current sheet and the strongly magnetized stripes separating the 
sheets - each sheet is modeled as a uniform unmagnetized slab lying on top of the 
confining uniform magnetic field in the effective gravity felt in the frame co-moving 
with the accelerating wind. Lyubarsky's chosen configuration is always unstable to 
incompressible "fluting" disturbances (in the sense of so-called flutes on a column), 
for sufficiently short wavelengths across B and long wavelength parallel to B. How- 
ever, as is well known ( |Arons & Lea| [l976, and references therein), the real inho- 
mogenous system, with magnetic pressure increasing from inside the plasma layers 
to the more rarefied and colder regions outside, is readily stabilized by the magnetic 
pressure gradient, against which potentially unstable "heavy" plasma mass elements 
have to do work as they try to fall into the low density magnetized region. Properly 
analyzing this instability requires a more sophisticated equilibrium model of the 
accelerating stripes and current sheets than that used in Lyubarsky (2010). 

Upon injection into the PWN, the flow decelerates at and beyond the wind "ter- 
mination working surface" (TWS) that occurs where the the dynamic pressure of 
the wind Er/4kRj WS c balances the pressure of the previously injected plasma and 
magnetic fields captured in the nebular bubbl^j In the Crab Nebula, the TWS is 
usually identified with the prominent ring of emission seen in the Chandra image of 
the inner nebula modeled as a MHD shock wave, located at about the same location 
as the other signs of activity seen as indicating the termination of the free wind from 
the pulsar ("wisps", etc., reviewed by Hester 2008 ) Because the lack of identifiable 
emission from the region interior to the TWS suggests the upstream flow is cold, and 
(according to MHD flow models for the relativistic wind) is super-magnetosonic, the 
TWS is thought to be the simplest flow termination "catastrophe", a relativistic mag- 
netosonic shock wave (e.g. |Kennel & CoronifI||1984[ ). Analogy to non-relativistic 
shocks in the supernova remnants, which clearly are efficient converters of flow en- 



ergy to non-thermal, power-law-like spectra of accelerated particles (e.g. Reynolds 



201 1| ), suggested to many that the non-thermal particle spectra inferred in the Crab 



and other PWNe arise from some form of shock acceleration. 

However, despite the sophistication achieved in MHD modeling of the inner 
Crab, especially of its variable wisp structures (Camu Fef al. | [2009 ), no clear identi- 
fication of the shock itself has appeared - in particular, the Chandra ring has eluded 
the MHD modelers. The applicability to other PWNe has also been unclear. At a 
basic physical level, it is unlikely that diffusive shock acceleration (DSA), assumed 



I discuss here only the young high voltage nebulae, where pressure due to motion through the 
interstellar medium and to reverse shocks from the surrounding supernova remnant do not have 
major effects on the structure - these young systems are the most useful for probing the particle 
acceleration physics and the plasma properties. 
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b) 




Fig. 10 a) Magnetic Geometry of a Force-Free Rotator for r < 2Rl, for i = 60°, from Spitkovsky 
(2006). The rapid transition to inclined split monopole field geometry for r > Rl is apparent, b) 
Geometry of the current sheet from the split monopole model for i = 60°, r> Rl. For clarity, only 
one of the two spirally wound current sheets is shown. As i — > 90°, the sheets almost completely 
enclose the star; for r^> Ri, the spirals are tightly wrapped (B r <C B§) and the current sheet surfaces 
closely approximate nested spheres, c) One sheet for ;' = 30°, shown for clarity, d) Meridional cross 
section of the current sheet for i = 60° . e) Equatorial cross section snapshot of the current sheet, 
showing the two arm spiral form. The arrows show the local directions of the magnetic field; 
the dots and crosses show the direction of the current flow. Panels b)-e) were constructed using 
Bogovalov's I Bogovalov 1999 ) analytic model of the asymptotic wind, f) Current sheet from a 2D 
PIC simulation of the inner wind, from unpublished work by Spitkovsky (used by permission), 



by most to be the microphysics in the shock acceleration, actually is at work at the 
TWS in the Crab and other PWNe. The magnetic field in the wind is transverse to the 
flow - in the Crab, the field winds up with ~ 10 9 turns between the light cylinder and 
the TWS. It has long been known that DSA has difficulties in transverse shock ge- 
ometry. In relativistic shocks, the process stalls because particles follow field lines, 
and as soon as the angle between the upstream field and the shock normal exceeds 
1 /Jl, with c/3ifi 3> c the upstream 4-velocity, the velocity normal to the shock of 
a particle attempting to return to the upstream medium becomes less than the up- 
stream fluid speed, so that bouncing between upstream and downstream becomes 
impossible (Begel man & Kirk"l|1990[ ) - the shocks are superluminal. That difficulty 
can be overcome, if large amplitude, short wavelength magnetic turbulence exists 
with amplitude sufficient to cause fast scattering across the average magnetic field 
( |Niemic & Ostrowski| [2006) |Lemoine et~al\ [20061 |Pelletier et al\ [2009). How- 
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ever, extensive PIC simulations of superluminal magnetosonic shocks in otherwise 
unstructured upstream media (no macroscopic upstream inhomogeneities, such as 
clumps in density, no current sheets, no particle component in the upstream flow 
with large Larmor radii at encounter with the shock,...) has shown that the shocks 
themselves generate no such turbulence ( |Sironi & Spitkovsky 2009 and references 
therein, Spitkovsky & Arons, in prepartstion), when the upstream magnetization 
<Ti = (B 2 /4»c 2 r) U p St ream is as large as is inferred in macroscopic MHD models 
of nebular surface brightness. 



2.4 Beyond MHD 

These and other problems have directed attention to extending the study of the TWS 
to conceptual realms larger than an MHD shock in an unstructured medium. The 
essential consideration, driven by the angular anisotropy manifest in the nebular 
emission around many pulsars, is that there is something special about the angular 
sector of the nebula around the plane of the rotational equator of the underlying star 
- the X-ray torus and the Chandra ring in the Crab and the partial X-ray tori in other 
PWNe ( |Ng & Romani||2004||Kargeltsev & Pavlov| [2008 ) indicate the significance 



of this region. The simplest special aspect follows from the angular dependence of 
the Poynting flux of a o" > 1 wind: c{E x B/An) ■ e r = (E R /4nr 2 c)cos 2 A, A = 
rotational latitude ( Michel} [l 973) >, a description that certainly covers the inner wind 
r <C Rtws- MHD models of the Crab Nebula's post TWS flow that use this energy 
flux as input yield good models of the X-ray surface brightness, if two not unrea- 
sonable additional assumptions are made: a) a particle flux uniform in latitude and 
b) an extra suppression of the energy flux toward the equatorial plane designed to 
represent hypothesized decay of the magnetic field upstream or at the TWS in the 
split monopole geometry of the outflow. For example, ( Del Zanna et al.\ 2006} con- 



sidered a total energy flux « cos 2 A tanh 2 (A /Ao), with Ao the angular half-thickness 
of the putative current layer across which the wind's toroidal field reverses. They 
found Ao ~ 6° and (oi) ~ 0.02 (the average is over An steradians) gives the best 
comparison of the synthetic surface brightness map to the Chandra X-ray images. 
Whether field decay occurs well upstream, as Coroniti envisioned, or occurs within 
the TWS itself (Lyubarsky 2003 ), remains an unresolved issue. Field decay within 
the TWS is of possible importance if the upstream magnetic stripes do survive all the 
way to the TWS, for then reconnection, driven by the shock induced increase in the 
pressure, can dissipate the oppositely directed stripes of magnetic field ( jLyubarsky] 
2003"1|Sironi & Spitkovsky] |20TT) 
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2.4.1 Current Starvation and Regeneration of Strong Vacuum-like Waves: 
Does it Happen? 

From a larger point of view, the non-trivial difficulties with the simplest variants 
of the standard model have, over the years, led various authors to suggest giving 
up MHD (even with dissipation) entirely. A hoary motivation to do so, going back 



to Usov (1975 i and (Michel 1994 1, is the claim that as the wind density drops 
oc r~ 2 , current in the stripes becomes starved, since /<x r 1 so that at large enough 
radius, typically estimated to be r > 100 AU s» 10 15 cm -C RjwSi J/2en± oc r ex- 
ceeds c, an impossibility for conduction currents. This idea comes from simple use 
of Ampere's law in the proper frame of the flow, with /' oc V' x B' « 1 /RjjT^^, 
from V' = V/I^md) V ~ 1/stripe wavelength = 1/Rl andS' =B/r w j m j. The (usu- 
ally unstated) assumption is that the relevant magnetic gradient scale is the stripe 



wavelength (e.g., Kundt & Krotscheck 1980 Melatos & Melrose 1996), with cur 



rent starvation setting in for R st arve < r <C Rtws- The supposition then has been 
advanced that MHD must fail, the frozen-in wave structure becomes a travelling 
wave ven in the proper frame of the plasma flow, and displacement current takes 
over from conduction current, thus resurrecting the strong, vacuum-like waves that 
were the foundation of pulsar rotational energy loss models in the early days of 



pulsar theory. Hypothesizing that this transition occurs rather abruptly, (Melatos 



1998) derived jump conditions describing the transition from the frozen-in, comov- 



ing standing wave (in the wind proper frame) described by the striped wind model 
(a nonlinear entropy wave, as perceived from the point of view of MHD wave the- 
ory) to the superluminally propagating, strong EM waves conjectured in this idea; 



following the same idea, Akra & Kirk (2012 1 rediscovered and extended these jump 



conditions, as well as substantially extending the equilibrium theory of strong waves 
in a cold plasma. 

Such models have not shown that the hypothesized mode conversion actually 
can be realized - that would require an evolutionary or instability calculation, show- 
ing that the nonlinear entropy mode actually can couple to the superluminal strong 
electromagnetic waves; the existing models assume that coupling occurs and use 
jump conditions to characterize the supposed transition. Viewed from the perspec- 
tive of small amplitude wave theory, that is a bit unlikely - the dispersion relation 
for the MHD entropy mode, co = (the / indicating plasma rest frame quantities 
is omitted for notational simplicity), never crosses that of the superluminal modes, 
except exactly at propagation cutoff, a non-propagating version of the wave where 



existing theory and simulation (e.g.,|Leboeuf et al. 1982, Skaeraasen et al. 2005 1 and 



references therein) suggests the waves are unstable and subject to strong damping. 
Nevertheless, the existence of even that limited resonance suggests that there might 
be merit in an evolutionary investigation, to see if the hypothesized transition might 
actually be realized - then a more properly posed stability investigation would yield 
more useful information on the plasma heating and radiation characteristics that 
could be useful in employing this kind of model for studying relativistic striped 
winds in the real world. A well chosen simulation is probably the best tool to study 
the question. 
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That idea that current starvation might occur and lead to a transition from MHD 
wind to some kind of superluminal electromagnetic wave was first made explicit 
by Usov ( 1975 i. His oft repeated argument (e.g. Michel 1994 ) may be ill founded, 



based as it is on the use of the stripe wavelength Rl as the scale for the current 
flow supporting the stripes. In fact, the stripes are supported by the current in the 
current sheet, whose thickness A probably is on the order of the formal Larmor 
radius r^ = m±c 2 (y±) jeB, where (y±) is the average Lorentz factor of the parti- 



cles inside a sheet and B is the magnetic field just outside the sheet (Michel 1994 1. 
Pressure equilibrium between the plasma in the current layer and the magnetic field 
outside yields ri = c/y/2co p ± - the current thickness is on the order of the skin 
depth, indicating the relation of this structure to an electromagnetic pulse penetrat- 
ing an overdense plasma. The magnetic field jumps across the current sheet from 
B(ji(= < I ) mag/ r ) to — B§. Integrating Ampere's law across the layer yields the surface 
current carried by the sheet j = cB§ /2n. Since j = JA, the volume current support- 
ing the stripes is J = cB^/lnA - the current sheet thickness A replaces the wave- 
length Rl that appeared in the conventional argument put forward by Usov and his 
successors. Assuming, as is plausible, that A is a factor k times rt, k ~ a few, I find 
J s» eB^2nkm±c(Y±) x r~ 2 - J/enc does not decline as the wind expands, except to 
the degree that plasma heating in the current sheet causes (y±) to increase with r. 
According to the this estimate, J/enc = B 2 ^ /2nkm±c 2 (y ± )n — B 2 ^ /2nkP± —4/k<l 
if the current sheet is at least 4 formal Larmor radii thick. If the MHD model satisfies 
this constraint, the stripes' fate is determined by magnetic dissipation, as Coroniti 
hypothesized. Whether that increase occurs requires properly accounting for the 
density in the current layer, which consumes the plasma in the stripes, and account- 
ing for radiative cooling as well as the finding a proper model of the dissipative 
heating. Most likely, if the wind was not current starved as it was launched, it never 
becomes so. But a proper determination awaits much more refined current layer 
models than have appeared so far. Since radiation losses offer the prospect of ac- 
tually observing something, such modeling should be of more than idle theoretical 
interest. 

Whether the dissipation does work fast enough to destroy the stripes interior to 
the TWS depends on physics not fully characterized so far. All that is certain is that 



if the mass flux is as small as that described by Kennel & Coroniti ( 1984 1, there 



isn't enough time to annihilate the stripes ( |Lyubarsky & Kirk|2001) > no matter what 
the detailed microphysics may be, but as discussed above, incorporation of the total 
mass loading, including that required to feed the radio nebulae (explicitly ignored 
by Kennel & Coroniti), does allow for reconnection velocities to be less than c, thus 
reopening the possible relevance of upstream dissipation. 



2.4.2 Observing the Wind 

This subject has an overly theological flavor, because of the lack of direct observa- 
tions - the theoretical physics is interesting, but as with all non-linear plasma prob- 
lems, making serious progress and filtering out the extraneous from the relevant 
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ideas needs experiments - numerical, laboratory and telescopic observations. The 
upstream wind and its putative short-wavelength structure cries out for simulation 
and, if possible, experiment, perhaps possible in laser generated plasmas. Even more 
important, detection and characterization of photons (or their lack, through absorp- 
tion of radiation coming from the pulsar) that can be unambiguously attributed to 
the wind would greatly advance the subject. Inverse Compton scattering due to the 
wind's supposedly relativistic bulk flow may have sufficient luminosity to allow de- 
tection of inverse Compton emission from the wind Bogovalov & Aharonian (2000 ) 
at r < 1007?^, perhaps distinguishable from brighter unpulsed nebular emission by 
its different spectral characteristics. 

The dissipation of the current sheets in the inner wind (r <C Rts, while respon- 
sible for only a small fraction of the magnetic destruction, might appear as pulsed 
emission to the observer at the pulsar's rotation period. The point is, if emission is 
confined to a sequence of thin layers separated radially by spacing R sep (= the light 
cylinder distance Ri in the nominal striped wind model), moving out with speed 
cfiwind = c ~ c /2r^i n( ii r w ind » T then so long as the emission at a given frequency 
has duration (in the pulsar frame) small compared to (RL/2c)r^ jnd = {P/47t)rl ml , 
the radiation from the sheets will not overlap and will arrive at the observer as a 
series of pulses occurring at the underlying pulsar's rotation frequency P. This al- 
ternate to the usual lighthouse beaming model for pulsed emission, first alluded to 
by Michel ( 1971} and more explicitly by | Arons ( 1979} , hasn't been explored much. 
The lighthouse model for pulse formation, was adopted early on because of its sim- 
ple explanation of the variation of pulsars' radio polarization with rotation phase, an 
explanation that has survived all challenges, then later adapted to gamma ray emis- 
sion. Gamma rays are now clearly understood to arise at radii much larger than those 
for the radio emission, which comes from deep within the stars' magnetospheres. 
The wind emission idea offers the promise of probing the otherwise unseen wind 
upstream of the TWS, so if such emission could be unambiguously identified, a lot 
could be be learned about the (experimentally) mysterious wind. The upsurge of in- 
terest in current sheets in the wind that followed Coroniti's model for the wind's B 
field dissipation has motivated some recent development of the idea of pulsed emis- 



sion from the sheets, by Kirk, Petri and collaborators ( Kirk et al. 2002 Petri & Kirk 



2005 Petri 2009 201 1 ) as an explanation for pulsed optical, X-ray and gamma ray 
emission from the Crab and other pulsars, while Aharonian et al. ( |2012 1 have of- 
fered a variantt of the Kirk et al\ ( |2002| l model with inverse Compton radiation 
emitted at r ~ 5QRl as a model for 100 GeV pulsed emission recently discovered 
from the Crab pulsar ( Aliu et al. 201 1| Aleksic et al. 201 1 1. That model should be 



compared to the inverse Compton scattering model with emission originating from 
near but within the light cylinder offered by |Lyutikov et al. (2012 1. 

So far all these interesting ideas have not found unique and incontrovertible use 
in modeling pulsar observation^]- there is huge theoretical freedom in this kind of 
phenomenological modeling, and so far there have not been any observations which 
uniquely require wind emission as the model (although the £ > 100 GeV pulses 



8 A criticism that also can be leveled at the lighthouse model for high energy pulsed emission ■ 
while more popular, it's not clear that it is is unique. 
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from the Crab come close to demanding consideration of the wind, as do the Bai 



& Spitkovsky (2010) lighthouse beaming models for the lower energy gamma ray 



pulses.) Parameterized model results of pulsed emission from the wind can be fit 
to observations of the light curves; however, so far they have had even less predic- 
tive predictive power than lighthouse models, which rely on relativistic flow along 
B rather than across B. Lighthouse models for outer magnetosphere emission do a 
little better for the pulsed optical emission, in the few cases when it has been de- 
tected, it that the magnetic geometry employed doesn't require ad hoc additions to 
the model, something that is required in the striped wind emission model - getting 
the sweep of optical polarization in the Crab required adding an ad hoc meridional 
component of B that rotates through the current sheet, a component not apparent in 
the force-free simulations of the inner wind. But current sheets could, in principle, 
have such embedded rotating vectors. Not having been designed to study current 
sheet structure, the existing numerical models may or may not offer a serious objec- 
tion to this interesting phenomenological model. More focused simulations would 
be useful. And in general, developing more self-consistent versions of these many 
interesting ideas would make comparisons to observations less subject to theorists' 
wiggling out of discrepancies by freely adding ad hoc components to the models. 

Detection of absorption and other radiation transfer phenomena in the light 
curves and spectra of the underlying pulsar, which backlights the wind, might offer 
a probe with less modeling freedom. There is a large literature exploring the influ- 
ence of propagation effects, including absorption, on the formation of the complex 
light curves of radio pulsars, both as alternates to, and resolutions of, difficulties and 
issues in modeling the frequency dependent light curves solely in terms of emission 
and geometric beaming characteristics (e.g. Wang, L ai & Han|2010| and references 
therein). There has been little such effort in using the properties of photon transfer 
through the relativistic wind in the interpretation of pulsars' appearance, at any pho- 
ton energy, although a few instances of such efforts can be found. For example, long 
ago, [Wilson & Rees] ( 1978] ) used the transparency of the wind to induced Compton 
scattering of the pulsed radio emission as a means of setting a lower limit on the 



wind's flow Lorentz factor in the inner wind they found rj, ,-, K / > 10 4 . |Barnard| (]1986 ) 
showed that the polarization limiting radius for radio waves of short period pulsars 
might fall in the inner wind, where the B field is toroidal and does not change its di- 
rection with pulsar rotation, thus plausibly explaining the lack of sweep of the Crab 
pulsar's radio polarization with rotation phase, even though the optical pulse does 
show rotational sweep of the polarization - the temporal coincidence of the pulses 
suggests they arise at the same place, where the B vector does rotate. These exam- 
ples show that transfer of radio waves through the wind might create imprints on 
the observed emission which could be useful in diagnosing the wind - radio waves 
interact much more strongly with the wind's plasma than is the case for the higher 
frequency optical, X- and 7 rays. 
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3 Flares: Analysis of Consequences 

The recent discovery of gamma ray "flares" in the Crab Nebula's emission (the data 
are discussed in £JTJ may have made the problems with the standard model more 
urgent. Not only is it hard to see how DSA can succeed in a shocked flow with mag- 
netic field transverse to the flow - it is hard to see how any accelerator with electric 
field less than the magnetic field, as is true of all models based on scattering of par- 
ticles in magnetic turbulence, can account for the large amplitudes, short time scales 
and high photon energies of the flaring events, either with or without an embedded 
shock. At photon energies above 100 MeV, the 0.1 — 1 PeV particles have very rapid 
radiative losses. The synchrotron cooling time for particles of energy E = ym±c 2 is 
(assuming a uniform distribution of pitch angles) 

6nm±c 2 9 c 28.4hours 

Tsynch = ca Thomson B 2 y = 4yi2^^ = B 3/2 1/2 / (12) 

muliGauss GeV 

Here Q c ± = eB/m±c 1 Oj] mmon = (87C/3)r 2 , r e = e 2 /m±c 2 , and the synchrotron 
radiation equivalence of the observed photon energy £ to the characteristic en- 
ergy e c = (3/2)hQ c ±Y 2 has been used to obtain y = ^(2/3)(s/hQ c ±) = 7.6 x 
10 9 £Gev/B m uiiOauss- Rapid scattering requires deflection of the particles in the 
magnetic field, which can't happen any faster than the relativistic cyclotron time 

T C yc = 2ny/Q c± = 754(E PeV /B mimGauss ) hours : (13) 

the GeV photon emitting particles lose their energy within a few percent of one 
cyclotron orbij^] Thus 

T cvc 87tn c± r e , 87T Q c ±r e (2/3)e I6n e 
T synch 9 c 9 c TiLl c± 27 m±c z 

independent of the magnetic field strength. (Xf is the fine structure constant. Radia- 
tion losses can occur in time short compared to the relativistic cyclotron time since 
emission of the photons in question occurs without the completion of a full gyration 
around the particle's guiding center - only a fraction ~ y 1 °f tne full Larmor orbit 
need be executed for a photon's emission. The time to radiate a photon is only the 
nonrelativistic cyclotron time 2n/Q c ±, not T cyc . 

Of course, in most astrophysical circumstances the radiation time is longer than 
the gyration time - usually, the particle has to emit lots of photons before its energy 
changes substantially. For the particular parameters of relevance to the Crab flares, 



expression ( 14 1 shows that the cooling time is longer than the gyration period for 



particles emitting at energies less than 30 MeV (E < 200 TeV). 



This estimate assumes the fluctuating magnetic field causing scattering and the overall magnetic 
field causing the accelerated, radiating motion (which includes the fluctuating field) have the same 
magnitude, as is assumed in models which assume diffusion at the Bohm rate. 
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The fact that particles emitting photons with energy above 100 MeV cannot com- 
plete full Larmor orbits makes traditional DSA, and also second order Fermi accel- 
eration, in the TS region unlikely as the basic accelerator at these highest energies 
- these schemes have particles bouncing between major deflections in magnetic tur- 
bulence, with the bouncing caused by their non radiative Larmor gyration in the 
up and downstream magnetic fluctuations. The bouncing is slow compared to the 
particles' Larmor times. Therefore, this mechanism cannot bring the particles up to 
the energies required for the observed y-ray emission. Independent of the magnetic 
field strength, accelerating electrons and positrons up to PeV energies requires an 
acceleration rate larger than the Larmor gyration rate. 

An alternate view of this issue comes from considering the maximum energy a 
particle can achieve in any accelerator in the face of synchrotron losses. The accel- 
erator has an electric field #, which may be laminar or stochastic, thus particles gain 
energy at the rate 

ace $ 

7 <? = n c± -, (15) 



and 



2 ""cat D : 
m±c A B 

v R T. R 

(16) 



y _ y B Tcy C B 



y Q c± S In S 
If radiation losses occur at the synchrotron rata^jwith loss time as in ( 12 1, then 



W = T cyc B = 4 r e Q c ± ^B = 8 e ^B 
T S ynch 2nT sy „ ch S 9 c £21 m±c 2 S 



As particles accelerate, the radiation loss time decreases, until at high enough en- 
ergy, the radiative losses balance the energy gains and no further acceleration is 
possible. Thus the radiated photon spectrum should show exponential roll-off above 



the characteristic energy e max obtained by setting the ratio in ( 17 1 equal to unity. 



27 m±c 2 B B 
emax = y g = 236- MeV, (18) 

an extension of the long familiar estimate (Guil bert it et al.| |1983| [de Jager et al.\ 
1996) to explicitly include the dependence on the accelerating electric field. Figures 



2]and[3]show that the flare durations exceed the synchrotron cooling time, therefore 
the accelerator is "on" during a flare - the radiation is not a consequence of an 
impulse of acceleration followed by cooling with S switched off during the emission 
- and that the particle energies have reached the radiation reaction limit. 

The spectra of the biggest flare to date have exponential cutoff energies somewhat 
too high to be interpretable as coming from an accelerator with § < B - the cutoff 
energy exceeds (fTHj) by a noticeable amount, a fact which has led to the suggestion 



( |Uzdensky et a/.||201 1 Cerutti et al. 2012 1 that the electric field is associated with 



magnetic reconnection in a region with magnetic topology similar to a sheet pinch 



10 That rate assumes only that the radius of curvature of a particle's orbit be the relativistic Larmor 
radius, not that particles actually complete full cyclotron orbits. 
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or a cylindrical pinch, where S > B in the middle of the pinch. This kind of con- 
figuration is known to lead to flaring acceleration in magnetic confinement devices 
(during "disruptive instability"), in the solar corona (solar flares), and in planetary 
magnetospheres (terrestrial magnetospheric substorms), where particles in the weak 
B region in the middle of the pinch can undergo runaway acceleration in the electric 
field rather than drift motions in the crossed S and B fields. However, this hint from 
the high cutoff energy might be mitigated if the source region undergoes relativistic 
bulk motion toward us during flare times - a Doppler boost factor ~ 2 would al- 
low the cutoff energy in the plasma rest frame to be consistent with S < B in that 
frame. Doppler boosts also relax the constraints of rapid time variability on the size 
of the region in question. The short time scales of the flares - both duration and 
even shorter internal variations, as seen in Figure [3] - may well be attributable to 
disruptive instability (such as tearing) of a current layer 

However, a Doppler boost of this magnitude is not apparent in the regions where 
high angular resolution observations at lower photon energies have revealed vari- 
ability. The "wisps", the "knot", the "anvil", the inner regions of the "jet" all show 
apparent motions with speeds less than 0.5c. Komissarov & Lyuti kov| ( |2011[ ) sug- 
gested that relativistic flow from the part of the curved TWS that is tangent to the 
line of line of sight, modeled as a relativistic MHD shock wave (a very nice ex- 



planation of the nebular "knot" seen next to the pulsar put forth by Komissarov & 



Lyubarsky 2003), could be a flare site with such bulk Doppler boosting properties. 



However, as modeled within the confines of MHD, the flow velocities aren't high 



enough to give enough Doppler boost. Also, to date, high resolution radio (Lobanov 



et al. 2011 1, optical and infrared observations (C. Max and R. Romani, personal 



communications) of the knot have shown no changes temporally coincident with 
the gamma ray flares. That lack of coincident change has been used to exclude the 
pulsar itself as the source of the flares; the same reasoning suggests the knot also is 
not the prime suspect. 

Bykov et al. ( 201 \) take the more conservative view that all we are seeing is in- 



terim ttency in the magnetized turbulence required in the traditional DSA model. The 
short time flaring time scales and the long time between major flares are attributed 
to assumed intermittency in the turbulence, with luminosity increasing as magnetic 
field increases and the hardening of the spectrum the result of the increased emis- 
sivity at E max causing the exponentially declining spectrum at higher energies to be 
visible at higher e than in quiescence - since radiation reaction limits the particle 
energies, the characteristic energy of exponential rolloff is independent of the field 
strength, therefore the shape of SED of the emission does not change during in- 
creases and decreases of B. Aside from the difficulty of making DSA function in the 
face of the strong radiation losses, the SEDs of the April 2011 flare look too hard 
to be readily explained by this idea - the increasing flux at 100-300 MeV seen in 
panels 7 and 8 of Figure [4] are not reproduced by the "toy" model spectra shown by 
Bykov et al. ( 201 1) . Nevertheless, variable B in the flare site surely contributes to 
the variability of the emission, on top of the impulsive behavior of the accelerator. 

My personal bet is that the flare site is the equatorial current sheet in the imme- 
diate upstream (r < Rjws) an d Chandra Ring region of the wind and TWS, where 
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the outflow velocity is still relativistic, and the current sheet separating the northern 
and southern hemispheres of the outflow is a permanent feature of the structural 
Tearing in the reasonably flat equatorial current layer destroys azimuthal symmetry, 
leaving an array of radial current beams and reconnected islands, themselves subject 
to kinking, causing Doppler beaming to have highly variable directions, a plausible 
source of the observed variability. The long build up to major tearing disruption, 
seen in the Earth's magnetotail (e.g. Ohtani et al. 1992), may be the origin of the 
long time between flares - that buildup time is limited by the transit time from pulsar 
to TWS, ~ months; minor tearing events, occurring quasi-continuously, may be the 
origin of the restless behavior of the light curve shown in Figure|2] Variable beaming 
also suggests we don't see all the flares that occur. The radial reconnection electric 
field has a large length over which it can do its work, allowing particles to reach 
energies up to a large fraction of the total wind voltage mag - this contrasts with 
the more limited acceleration possible in the tearing of the closely spaced current 
layers in the equatorial flow region, when the stripes are assumed to survive all the 
to the TWS, as in the model studied by Sironi & Spitkovsky (2011). The config- 
uration is that of a linear accelerator (as in the earlier work on a similar, AGN jet 
motivated model by Larrabee, Lov elace & Romanova|2003 and references therein), 
therefore radiation reaction limitations are reduced, especially for particles which 



stay well focussed within the current layer (Uzdensky et al. 2011 Cerutti et al. 



2012} - the current layer is relatively thick, with characteristic half width set by the 



latitudinal toroidal field gradient set by the imbalance in the magnetic flux in the 
neighboring stripes at latitudes off the exact mid-plane (Coroniti"} |1990 1, thus longer 
particle residence in the current layer and higher voltage drops are possible than is 
the case where the current sheet forms oppositely directed current layers situated 
only Rl apart from each other, as in the striped geometr^p^] A configuration without 
the short wavelength = Rl structure, such as the equatorial current layer remaining 
after upstream decay of the stripes, has a huge advantage as a geometry for recon- 
nection powered runaway acceleration in the current layer. The accelerated particles 
form beam dumps at the ends of the current filaments, perhaps providing a model 
for the hot spots observed on the Chandra ring. Because the particles are runaways, 
the particles form a beam - in the limit of a strictly ID beam model, the particle 



Coroniti 



1990 



The 



11 The picture I favor is close to Coroniti's model of dissipation in the wind ( 
striped magnetic structure launched from the pulsar (see Figure \W) decays interior to the TWS, 
probably at r ~ 0. IRjws (Arons, in preparation), which is possible in the highly mass loaded wind 
in the Crab and other young pulsars. 

12 The focussing of particles toward the current sheet's center is a rediscovery of the focussing 
principle long known to accelerator physicists (e.g., Courant & Snyder 1958 ). That mechanism can 
be used as a means to confine particles over distances ~ 10' 'Rl = sheet spacing, a length required 
if voltages of a PetaVolt or more are to be accessed by the linear accelerator, only if the fields 
are smooth and properly "designed" to a degree extraordinary for an astrophysical configuration, 
especially when subject to the uncontrolled macroscopic gradients introduced by the instabilities of 
the current sheet. The Sironi & Spitkovsky simulations of the striped configuration show that when 
the formal Larmor radius of the particles in the sheets becomes somewhat larger than the spacing, 
the acceleration saturates, because particles have finite Larmor radius drifts out of the sheet's core 
in the "messy" fields of the unstable layer. 
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distribution is mono-energetic, but in general, the spectrum should be very hard, 
perhaps even inverted in energy space, thus explaining both the appearance of a new 
very hard component in the gamma ray spectrum during a flare, and the negligible 
emission at photon energies below 100 MeV). This scenario suggests the emission 
region to be the outer wind, at angular separation from the pulsar ~ 0.2" — 5", but 
with faint emission at energies below the VHE band that twinkles, because of vari- 
able Doppler beaming, That time variability suggests the value of continuous high 
sensitivity and high angular resolution monitoring in radio, optical/infrared and X- 
ray bands - arc second angular resolution in the VHE band is out of the question, 
unfortunately. A quantification of these ideas will emerge soon. 
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